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PREFACE 


This paper is intended to provide the theoretical and practical 
framework for a two-dimensional temperature model. The motivation i 
for this study is provided by the existence, at Rand, of a general 
two-dimensional flow-reaction simulation modei ueveloped by J. 

Leendertse and £. Gritton. The vertically integrated mass balance 
equations of the simulation model are discussed and shown to apply 
to temperature modeling. Terms contributing tu the total heat flux 
across che water surface are described and approximated by empirical 
relations wnich, in turn, are integrated into the temperature model 
through an equilibrium temperature and heat exchange coefficient. 
The sensitivity of these parameters to the meteorological variables 
is discussed in the context of a preliminary calculation. Tables 
are given for the heat exchange coefficient and equilibrium tem- 
perature. Finally, expressions are derived for the relative iso- 
baric slope as a function of chlorinity and temperature. 

Although this paper is primarily a review, the approach of 
Edinger and Geyer is modified for the case when data is averaged 
over periods of less than a day. The following advantages of the 
E-G approach appear to remain: 


(1) Additional arrays to compute the surface heat 
exchange at each grid point are not required (although 
it may be necessary to divide the bay into large sub- 
regions with different Tt and K). 


(2) The method is readily adaptable to tong term cal- 
culations. 


Sections 3 and 4 of this report (on computation of absorbed 
radiation) were in. luded to emphasize the need for further research. 
The empirical formulatifons discussed are essentially useless in 
the development of a temperature model, but improved verions of 


these equations could be quite valuable in tapping the model's 


Te Tn eeaahnenanal 


predictive capacities. 
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IT. THE HEAT TRANSFER EQUATION 


It is desired to exploit as far as possible tne similarity between 


the heat and mass transfer equations. The former may be written: 


ot + Os (uT) i. 3 (vT) “gs a9. (WT) | Oo K ot 
ot ox dy 3z 3x x Ux 

2 nas 7 y S (1.1) 
pig EU), geet 2) pees ; 

ay y jy az zZ4Zz ia 


where p = density of water, assumed constant at 62.4 ib-ft°. 


hehehe = ot beatae eememen C . specific heat of water taken 
as 4.1 Joul-g ~ °C ~ = .98 Btu-lb ~ °F (assumed constant although it 
may vary as much as 17 due to salinity and temperature distributions 
in the bays under consideration, and far more in some estuaries), So = 
local heat source and sink. 

Next, introduce the following distributions and vertically aver- 


aged variables: 


t 

ae Ll ue 1 

T = A ff ta: eT €T)° 5 (1.2) 
-h 

T(z) = Tlie. (z)) , (1.3) 

u = U {1.- cy (z)] , ete. (1.4) 


(1) 


The exact vertically integrated continuity equation is 


23 #9, Cy = 0, where v= (u,v) and e = 
(1.5) 
ae 
ax dy] 
and the boundary conditions on the free surface and on the bottom are 


wo ces v > Fr > Ww Eo Giang mee eh. (1.6) 


z jt z ch 9t hy 


Vertically integrating (1.1) and substituting equations (1.5 - 


1.6), we obtain 


3 oT 
/oc — 
Qy/ee, + («, 2) + (ay ( ‘ 3) CE x aE =) 


- (2, - 7, %- T oht y [ #. Cet) 2a, Wo = T. nl ni] 


e3t sh Oey 
A) > 
- - : = — v 
Pate Wt ac (T) + ¥. 7) (1.7) 
Here, Q, is the net heat flux into the vertical water column 
through the surface and bottom. In terms of the quantities Oe and 
(2) 


ar first introduced by Leendertse, where 


Gg = + (1+ e, (2) ep (2)) , ete. , (1.8) 


the right hand side becomes 


a. (iT) 


¥ T V) (1.9) 


RHS +3 (a T U4 25 (H « 


uT ay vT 

If T is uniformly distributed over the .2rtical, then a = l. 
This assumption is adopted throughout this paper, but the general 
equation may be the subject of future investigations. The vertically 
integrated diffusion terms are discussed by Leendertse for the case 
of the mass balance equation. The same analysis is assumed to apply 


here. 


II. HEAT EXCHANGE COEFFICIENTS AND EQUILIBRIUM TEMPERATURE 


The Leendertse-Gritton general mass-balance model for n reacting 


constituents is described by the matrix equation 


3(HP) , 2(HUP) , a(HvP) _ a. (HD, 3P ) 
cf 3x oy oat OX 
= (2.1) 
3 (HD, 3P ) a ra 
ay y By + {K] HP + HS = QO 


where [K] = reaction matrix, and S = source or sink vector. 


For DO-BOD reactions, the K-matrix and § are written in the form 


Kos JH 0 
(2.2) 


{K} 
21 “2 


| 
ut 


(2.3) 
0 


This suggests an analogy with the linear heat exchange model of 


and the equilibrium tem- 


Edinger and Geyer, where K/pC_ replaces K 
7 (4-7) 


1 
perature replaces ers Here K is the bee exchange coefficient. 
Ky has no counterpart in the :-G model but can be utilized to describe 
the temperature dependence of p.iucesses such as reaeration. For com- 
putational and theoretical convenience, therefore, we shall employ the 


following expression for heat transported across the surface: 


Q = OT - T)K (2.4) 


The remainder of this section will be devoted to deriving expressions 
for K and ae in terms of meteorological variables. 


The energy budget equation is 


Gc Oe Qe 1G a OD (2:5) 


net anet 


Py ie 
where 

Q = net solar radiation 

net 
= atmospheric radiation 

Qanet P 

Qr = long-wave back radiation 

Q_ = convection heat exchange 

Q. = rate of heat transport due to evaporation 


The details of the energy budg’t terms will be discussed in later 
sections. Our concern here is with the dependence on water surface 


+ Q 
net anet 
To» the other terms may be written as follows: 


temperature. In particular, while (Q ) is independent of 


= 4 
Vr? Zo +T_) 
Q. = 0 (Ly - Ly T.) (a. - ef (2.6) 
Q. = 0 B(Ly = LT.) (TS - TE 


a 
vapor pressure in inches Hg, respectively; 


where e , e_ are the saturation vapor pressure at Ts and the ambient 
s 


p = water density, D = .1661 x ig” Basher ee 
Ty = 460°F 

Lg = 1093.2 Bru/Ib. 

L, = +57 Btu/1b.-°F , 


and f, B are meteorological variables. 

Whereas (2.4) implies that, without advection, T instantaneously 
decays toward UG with decay time Cons it is clear from (2.6) that 
considerable complexities are buried in K and Teg Specifically, *% 
is not independent of TS) and Tag is the solution of a transcendental 


equation. 
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Furthermore, T, appears in (2.6) rather than the vertically aver- 
aged temperature, T. Even in supposedly well-mixed cases, it may not 
be a good approximation to replace T, with T. Measurements of the 


vertical temperature distribution may be needed to settle this ques- 
‘ ba 6 
tion (and the related matter of determining aT and a 5s ) 


Substituting equations (2.6) into (2.4) we obtain 


4 4 
Qe = Qae » Qeaet? ae (T) . i) 
- 0 (Ly - Ly T)) (e. - ef -oB (Lo - Ly T)) (T - TF 
= K <e - T) (2.7) 


We shall assume that T. = T. Then eG is defined by the equation 


(Q + Q 


4 
net anet? 72 ts a Ted? 


- (Lo - L, Pn [e.g - e,) + oe -t)| of 


= Gtr.) = 0 (2.8) 


which says that ue is the surface temperature (= vertically averaged 
temperature) at which the heat transported across the surface is zero. 
Edinger and Geyer do not solve (2.8) directly. Rather, they 
first linearize the back scattering term and then assume the relation 
Ss - e,) / a - Ty) = (e. - e) / ‘T. - Ty): The latter relation 
is not a bad approximation if temperatures are averaged over a day 
or more, so that Ts anv Tr are nearlv equal. For our purposes, it 
is a very poor approximation which introduces an unwanted Ts dependence 
into ae and thereby necessitates computation of ee at each grid 
point. Our approach is to find a numerically by ‘Newton's method. 
This is quite effective since Cage is a monotonically decreasing 
function which can be expressed analvtically if the saturation vapor 
pressure at eg is evaluated with the Lamoreaux exponential approxi- 


mation. One may obtain the initial value of Ta for the iteration 


8) 
from the previous time step. The Lamoreaux approximation is given by ~ 
6 7482.6 
T) = 412 2 oS : 
e(T) 6.4129 x 10° exp ( 398.362 5 pD (2.9) 


| 
| 
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which, when combined with equation (2.8), yields the iteration formula: 


ToS, ee /G°(T,) (2.10a) 


a 
ates | 


where GT) is given by equation (2.8) and 


. 2 ky 3 . 7 
G°(1,,) = ~4D(Ty + T,.)? + Lye [(e,, - e,) + “ ,)| 


2.1 
se eit.) pave. =e es 
0 leq eq| (398.362 +T_ T. os | 
The procedure outlined above is quite efficient. If the equili- 
brium temperature is computed hourly so that it changes only a couple 
of degrees at each time step, no more than two iterations are required 
to update T_. 
eq 
Since the equilibrium temperature will characteristically pass 
through the water temperature, T, twice daily on its roughly sinusoidal 
path, it is evident that equation (2.7) is not well suited for direct 
calculation of K (i.e. the coefficient of K might be zero). 


Alternatively, one may subtract (2.8) from (2.7) setting Lo - ut. = 
constant = 1050, to obtain: 


4 4 
(cr, + Taq? - (Tp + T) + 1050 of [og -e)+ BG - T) 
= K(T,, - 1) =k (cr, +T)) - My + 1)| (2.11) 
or, 


3 2 
K =D [«t, Pi Oy Pi aye Bye et 


2 3 { 
(Tp + T)” + (Tp + T) | + 1050 pf ,R+ [eng - e,) / 


| 
(To - 1 | (2.12) 
eq { 


2 ‘= h - 
The term ee e,) / (T, T) may cause troubie when ee T) 


is near zero, but for ee - T}] <10°F it will be sufficiently accurate 


to use the relation 


aeeetneemiind Mt eats 


. ; Bere) piminiipanstaies sentinel 


3 late meine lp a 


me 
- Tels 
ie = e7) / oe -T) = oT me = 
2 /1+T 
173.004 = 
796.724 +T + i ‘ 2 (2.13) 


Since K depends on T, it would seem necessary to calculate it at 
each grid point. In section VII we see that K is not very sensitive 
to T, so that usually one value of K will suffice for the whole bay. 

The linear approximation to equation (2.11) devised by Edinger 


and Geyer amounts to replacing the first bracket by 4NT i and the 


0 
factor ae - e,) / oon - T) by (e. - e,) / (T, - Tp: Neither 
assumption is very good for T averaged over less than a day. 

If our working hypothesis, that T=T o> turns out to be untenable 
we need merely replace T by T, in (2.7), (2.11-2.12), and employ meas- 
ured values of TS: It may also be possible to set TS = T{lte,(o)}, 
where eq (2) is the vertical temperature distribution. 

Before describing the general features of ae and K, including 
sensitivities to fluctuations in the various terms of the energy 
budget, the next few sections will be devoted to calculating contri- 
butions to the total surface heat transport. 

We note in passing that an alternative scheme to that of Edinger 
and Geyer would be to use equation (2.7) for the heat flux Q, and in- 
sert it directly into equation (1.7). This would require computation 
of Q, at each grid point. The evaluation of K and ta is more complex 
than that of Qps but the former does not require new spatial arravs. A 
further advantage of the E-G scheme is that average values of K and 
Tg obtained from short term representative periods (e.g. tidal cycle) 
can be used in long term calculations since these values characterize 
the entire waterbody, whervas Q., is a local variable. In addition, 
possesses natural periodicity which is amenable to Fourier re- 


(3) 


presentation. 


a a 
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III. NET INSOLATION 


The solar radiation on a horizontal surface at sea level is pri- 
marily a function of solar altitude, cloud cover and cloud height. 
Tne solar altitude can be obtained from the cate, the hour and the 


latitude. In particular, we have 
sin o = sin ¢ sin 6 + cos 6 ccs $ cos B (3.1) 


where > = latitude, 46 = solar declination, and 8 = local hour angle 


{LHA). The local hour angle is given in radians by 


wT 
t 


= GG (Greenwich hour angle - longitude) 


a (Local standard time + No. Time Zones to Greenwich 


longitude/15 - A - 12) (3.2) 


Approximate values for A and § are tabulated in Table i, 62249) 


The total incident solar radiation may be obtained from the fol- 


(11) 


lowing empirical formula based on Raphael's curves: 


\ 
329 sin | [1.05076 eC RIOLG bs 115" uy? : : 


Q. = 
when sin « > 0; 
Qs = 0, when sin u < 03 
. = O, when vin a > 0.26; 
gos 702264, when’ O-< Sin > 0226. (333) 
q 


Oe as given by (3.3) gives sufficiently accurate values for clear 


skies. For cloudy skies these values must be modified by an appro- 


priate factor. The cloud factor given by Raphael and based on the 
work of Fritz is (1 - .0071 c*), where C is expressed in tenths of 

skv covered. Comparison with available data reveals thit this factor, 
when corbined with Me above, produces accurate results for mean 
monthly solar Padtatton: The negative bias observed by Anderson 


for fosby's factor does not appear here. However, as pointed out 


Feb 
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May 


Jun 


Jul 
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Dec 
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by Anderson, large deviations may occur for daily solar radiation 
unless this factor is further refined to include variations due to 
cloud height. It is essential to accumulate an extensive record 
comprising hourly pyranometer measurements and cloud cover and cloud 
neight readings. Althougo many workers would prefer to utilize 
field data rather than to concern themselves with the intracacies 

of cloud attenuation, the latter option may lead to greater flexi- 
bility and predictive capability due to the existence of long term 
meteorological records suitable for statistical analysis. The fol- 
lowing empirical formula for insolation on a cloudy day should be 


regarded as tentative: 
2 
Q, = (Gl - 0.0071 Cc’) Q, (3.4) 


The reflectivity of the water surface, R, is defined by 


Ree 0 tig Si (278) 


and 


Q = Q, (1 - R) (3.5b) 


where Cee is tne net solar radiation entering the water after re- 
flection at the surface. R is principally a function of solar alti- 
tude, type and amount of cloudiness, and atmospheric turbidity. 
Turbidity has been shown to be relatively unimportant except in 
notably extreme cases. Anderson has given the following general 


empirical formula for reflectivity: 


Ree ears (3.6) 


where a is solar altitude, as before, and a and b are experimentally 


determined parameters resated to tne amount and height of clouds. 


If the cloud dependency of a and b is ignored, errors of order 57 
or less are incurred in the total daily Q5- To avoid accumulation 
of error in computations for periods of a week or longer, the fol- 


lowing empirical formula, based on Anderson's curves, may be used: 


-b = 1.1073 Log, )(8/D) (3.7a) 


1.9031b 


1/2 (N/10° + D/10 ) (3. 7b) 


f 
Li} 


where 


we 
" 


(.2004 F + .2341 F + .2006 F, + .1172 F ) (3.7c) 
c s b fo) 


oO 
ti} 


(.04040 ie + .03115 Fy + .04877 F, + .04694 FY) (3.7d) 


b 


and Poe Poo Fi Fy are the frequencies respectively for clear skies 
(C=0), scattered clouds (C=1-5), broken clouds (C=6-9), and overcast 
(C=10). To avoid infinite reflectivities, set R equal to its value 
at 5° when 0 <s a < 5°. 

In summary, the net solar radiation entering the water after 
reflection at the surface is given by equation (3.4-3.46), i.e. 


2 . 
Qrer = (1 = 0.0071 ¢*) (1 - aes Q. (3.8) 


where a and b are given by equation (3.7) and Q, by equation (3.3) 


f= 


IV. LONG-WAVE RADIATION 


Two ways for obtaining the net atmospheric radiation will be 
discussed. One is direct measurement. The other is through an 
empirical formulation promulgated by Audevson: (-? 

Atmospneric radiation can be measured by subtracting pyranometer 
(shortwave) readings from pyradiometer (long and short wave) readings. 
As pointed out below, good measurements are prudent before acceptance 
of any empirical scheme. 

Not all factors determining the atmospheric radiation are known 
or quantitatively understood. The primary ones, undoubtedly, are 
the total moisture content of the atmosphere and the temperature of 
the constituent water vapor. Any attempt to by-pass these factors 
by using local vapor pressures and ambient temperatures is doomed to 
partial failure. In particular, Anderson's study revealed that a 5% 
to 10% error is inherent in such efforts. 

Anderson assumes that the atmosphere acts like a "gray body" of 
variable "grayness", obeying a modified Stefan-Boltzmann law: 


4 
Q, = Boa (Tp + TD) (4.1) 


where 5 = Stefan-Boltcwanr. constant = .1712x19 eBtu hr tet R 
and ; = atmospheric moisture factor = a+b ev (4.2) 


The parameters a and b are given by the following formulas: 


a -740 + .025 C exp(-.0584h)} (4. 3a) 


b 


.164 - .018 C exp(-.060h) (4. 3b) 


where C is the cloud cover in tenths of sky covered, and h is the 
cloud height 1n thousands of feet (1.6:° h-~- =). 

The partial success of Anderson's anproach hinges on the appro- 
ximate correlation between local vapor pressure end total moisture 
content of the atmosphere. Unfortunately, any parameter fit re- 
tlecting tnis relationsnip must be restricted to use in locations 
which have similar air mass characteristics. Failure to aetermine 


parameters locally may, tnerefore, introduce bias of unpredictable 
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Magnitude, in addition to the candom error of order 10% already 
present. 

In summary, Anderson's equations (4.1 - 4.2) may be used with 
reasonable accuracy provided they are calibrated by adjusting the 
nunerical constants in a and b for each new location. 

The long-wave radiation reflected from the water surface is 
approximately 3% of the impinging atmospheric radiation for the 


normal range of water temperatures. Therefore we may write 


Q POOP Oss © (4.4) 


anet 


The water body itself radiates like a gray body of emissivity 
0.97 provided no film (e.g. oil slick) is present. This value is 


essentially independent of salinity. Therefore, 


ra 4 
oie = 0.97 o (T + T) (4.5) 


so that D = (.97) (.1712x107°) = .1661x10°8 Btu nro teem2 R74. 


. a ee 


pe ee TE I wae 
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V. EVAPORATION 


The results of evaporation studies to date are not éericluntye, 0ro 


Anderson's exhaustive Lake Hefner investigation produced an empirical 


equation of the form 
Get EE: AED OE! Rew Ee? he (5.1) 
e rey lis s a 
where 


f=MW ft hr incie (5.2) 


and W = wind speed (m.p.h.) at a specified height, Ho» above the water 
surface; M = empirical constant; e, = saturation vapor pressure (in- 
Hg) at the water surface temperature; and e ambient vapor pressure 
of unmodified air (i.e. unmodified by passage over the waterbody). 


Using equation (2.9), e, is given by 


e.* e(T.) (5.3) 
while e. is obtained either from 


e. > e(T,) (5.4) 


where Ty is the dew point, or from 


e =e(T ,) - P(T -T.,) (3.59x10 4 + 2.34x10 ” T,) (5.5) 
a wb a wb w 


where Ta, = wet bulb temperature, aie = air temperature or dry bulb 
temperature, and P = atmospheric pressure in-Hg. 

Anderson found that M = 2.08x10"" ft hr? iniig In addition, 
the Lake Hefner work appeared to support the theoretical ettorts 8) 
of Sverdrup and Sutton. Later studies contradicted these theories, 
however, and there presently exists no reliable formula derived from 
either continuous or discontinuous mixing concepts. All that has 
survived the series of investigations conducted by the USCS is an 
empirical correlation between M and the size of the waterbody for 
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the case when Hy = 2 meters: 


215= 


= gti! fone ean ie? : (5.6) 
x -05 
{[Note: to convert M in (in day mn) to Min (ft ee fraigo) mul- 
tiply by .11759.] Here, A is the area of the waterbody in acres. The 
standard error of (5.6) is given as 16%. 

The reasons for the above correlation are clear in a qualitative 
sense. First, there is a downwind decrease in evaporation rates. 
Generally, this will decrease M for increasing area, and, in addition, 
introduce a fetch dependence into “I. 

Secondly, the wind profiles depend on the size of the waterbody. 
With increasing fetch, the wind shear resulting from surface drag 
decreases. For cxumple, if the wind speed at 8 meters were the same 
for a large and small lake, then at 2 meters the wind speed over the 
large lake would be greater. As a result, if one wishes to use wind 
measurements made near the water surface (e.g. at 2 meters) one must 
compensate by using a smaller M for larger water surface areas. 

Since (5.6) gives only a fair approximation to these effects even 
for regularly shaped lakes, this equation should be applied to geomet- 
rically complex estuaries with some trepidation. It is quite likely 
that for Tampa Bay evaporation rates depend in a complicated fashion 
on location in the bay and on wind direction, as well as on wind speed. 

Two final points should be made regarding the use of equation 
(5.6). This equation applies only if the anemometer readings are 
taken in mid-bay. Also, it applies only to measurements at 2 meters, 
a height chosen to minimize the effects of atmospheric instability. 

If some measurements are made at different heights, the wind speed 
input must be adjusted accordingly. For example, the wind speed at 
4 meters is approximately 67% greater than at 2 meters for a 50,000 
acre bay. If a wide range of heights are used, detailed knowledge of 
the wind profile mav be required to obtain tne appropriate correction 


factors. 
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T. CONVECTIVE HEAT EXCHANGE 


There is no easy way to measure the convective heat transport nor 


any reliable way to calculate it. Consequently, one must rely, how- 


ever reluctantly, on a 45-year-old prescription due to Boweny 202)? 928) 
Q. = B(T, - TS) / (e, - e.) Q. (6.1) 
where 
-4 op 7l 
B = 3.39x10 P in. Hg.°F (6.2) 


and P = atmospheric pressure in inches Hg. 

The physical content of (6.1) is essentially that the evaporation 
eddy diffusivity ~ °~ the eddy conductivity are equal. 

Generally, one can ignore conduction in comparison with convection 
in much the same fashion as one ignores the contribution of molecular 
diffusion to evaporation (perhaps unless a strong thermal gradient 
increases free convection to a significant degree). One further assumes 
that the effects of spray on evaporation and the effects of radiative 
transfer (in moist air above the water) on conduction are negligible. 

The equality of eddy coefficients has been questioned for some 
time and several researchers nave modified the Bowen constant in in- 
stances where unstable atmospheric conditions prevail. The magnitude 
of deviation is thoroughly speculative, however, especially for large 


(5) 


lapse rates associated with high temperature outfalls. ‘ Even the 
careful work of Anderson did not reveal a patterned response for B to 
unstable conditions. Consequently, we shall adopt (6.1) and (6.2) 


without change. 


<7 


VII. PRELIMINARY CALCULATIONS 


In order to get a feeling for tne equilibrium temperature and heat 


, 


exchange coefficient, we shall calculate ne and K for a typical August 
day in Tampa, Florida. A useful heuristic tool for discussing the 
(4) 


response to se and K is the one-dimensicnal model o: Duttweiler. 


He considers the non-advective heat equation 


3T 
LE eee = (7.1 
v a T) rf ) 
with toe expressed in a Fourier series and k a constant. Taxing 
= + i + ° 
Tak? T T, sin(2nvt + $) (7.2) 
where T_ = time average value of Tas T, = amplitude of variation of 


Teg v = frequency, and ¢ = phase angle, he obtains 


2| -1/2 


T(t) = Tt T, [, + (27vh/k) sin(2mvt + 9 - 2) (7.3) 


+ F(O) exp(-kt/h) 


where v = tan} (2nvh/k), O<a<a1/2, and F(O) is a transient depending on 
initial conditions. Assuming that the model makes sense qualitatively, 
we find for a period of 24 hours (= 1/v) with k = a0." ft/sec and 
h=10ft that the amplitude variation of the equilibrium temperature is 
about 25 times the amplitude variation of the water temperature. Also, 
since a=m/2 for short periods like a day, one finds that the maximum 
water temperature lags a/2nv = 6 hrs. behind the maximum equilibrium 
temperature. 

An important consequence of these remarks is that the daily period- 
icity of T can be ignored insofar as K is concerned. This follows be- 
cause K is essentially a function of (T + Fa and the variation of T 
is negligible compared to that of Tae Of course, we must also consider 


the effect of spatial differences in T, particularly in the region near 


a high temperature outfall where K might be higher than in the rest of 
tne bay. 
Before dealing with this question, we calculate K and a from 


the following data: 
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Average Meteorological Input 
Solar declination = 14°N 

Maximum solar altitude = 76° 
Cloud cover = 6 tenths (C=6) 
Cloud height = 7800 ft (h=7.8) 
Cloud cover frequencies: F 


c 


F 
° 


1.3%, Fs = 23%, FL = 51.1%, 


UT 


24.4% 
Wind speed, W = 9m.p.h. 

Max (min) air temperature, T. = 90/74° 
Vapor pressure = .819 in Hg 


Dew point temperature, Ty = 73°F 


Bowen constant, B = .0102 in Hg opt 
Water temperature, T = 83°F 


From equations (3.8) and (4.4) one obtains the maximum and minimum 


radiation flux: 


(Quer + Wnep) = 365-5/123 Btu a a a eye oe (7.4) 


net a 


Substituting the radiation flux and the meteorological variables into 
equations (2.8) and (2.10), we obtain Famer) = 103.6°F after three 
iterations (beginning with a default value of 120°F) and tnt) = 
71.4°F after two iterations (beginning with 70°F). Therefore, the max 
and min equilibrium temperatures for the entire bay are 103.6 and 71.4 
respectively. From equation (2.12), the corresponding max and min 


values of the heat coefficient K are 9.23 and 6.68 Btu fice aes op) 


Tnese yield values .419x10- and .303x107" ft/sec for Ray = K/pC_. 

Nearly all the difference in K between the maximum and alain 
temperatures is due to the increased evaporation associated with the 
increased vapor pressure at the higher temperature. Verv little change i 
in K is caused by augmented back radiation. For the tvpical meteoro- 


logical conditions considered, the percentage variation in K per degree 


temperature variation is approximately “4K/(K4T) = 1”/ °F. This value 


~-j19- 


declines to a minimum of 1/3% op l when evaporation is zero, at which 
point it is due to back radiation alone. 

The previous remarks apply to variation of either Me or T. It 
appears therefore that high temperature outfalls wovld generate sub- 
Stantial deviation (10-20%) in K. Fortunacely, most of this deviation 
is of no consequence. In particular, as stated by Edinger and Pritchard, 
surface cooling is of relatively little importance in the initial and 
intermediate mixing zones. Only a small fraction of the discharged 
heat is dissipated into the atmosphere within the region where the 
temperature excess is 20% or more of the initial temperature excess. 
For example, if the water near the outfall is 20°F higher than the 
average for the waterbody (the temperature difference is usually 
smaller), most cooling occurs in water which has been diluted to less 
than 4° above average. Thus, in the area of principal concern K is 
less than 4% above its average for the bay. This discrepancy is 
probably less than errors involved in evaluating evaporation and 
conduction. The percentage variation in K per m.p.h. wind speed 
variation is AK/KAW = 100/W for typical evaporation rates. Therefore, 
if W is measured at 10 m.p.h. with 1/2 m.p.h. precision, then K is 
known with 5% accuracy. Values are given for the heat exchange co- 
efficient and equilibrium temperature in tables 2 and 3 respectively. 
Note that equilibrium temperatures above the boiling point of water 
occur frequently for W=0. This may be taken to signify that molecular 
diffusion becomes important in those instances. One should therefore 
set W equal to some small value (e.g. take W = .4m.p.h.) to 


min 
account effectively for free convection. 


(6) 


Heat exchange coefficient in Btu per hour per square foot 

per degree Fahrenheit as a function of water temperature 

=T, equilibrium temperature =E, and wind speed =W, 

Mass transfer coefficient is M = .000231 corresponding 

to wind measurements at 2 meters in waterbody of 50,000 acres. 


WATER TEMPERATURE = 55.0 DEGREES FAHRENHEIT 


W=0 W=2.5 W=5.0 W=7.5 W=10 W=15 W=20 W=25  W=30 
Eq Temp 


40.9 287 1.72 2.58 3.45% 4.29 6.00 7.72 9.43 11.14 
45.0 88 1.78 %.67 3.57 4.47 6,26 8.05 9.85 11.64 
50.0 -89 1.83 2.77 3.71 4.65 6,53 8.41 10.29 12.17 
55.0 91 1,89 2.88 3.87 4.85 6,82 8.80 10.77 12,74 
60.0 292 1.96 2.99 §,03  $.06 7.1% 9.21 11,28 13,35 
65.0 093 =62,02 3.11 4.20 5229 7.47 9,65 11,82 14,00 
70.0 095 2.09 3.24 4,38 5.53 7.82 10.11 12,40 14,70 
75.0 296 2.17 3.37 4,58 3279 8.20 10.61 13,02 15,44 
80.0 098 2,26 3.55 4,84 6.13 8.71 11.29 13,86 16.44 
85.0 99 2.36 3.72 5.09 6.46 9,19 11.93 14.66 17,39 i 
90.0 1,00 2.46 3.91 336 6.81 9.72 12.62 15,53 18,43 : 
95.0 1.02 2.56 4.11 5.65 7.20 10.29 13,38 16.47 19.56 
100.9 1.03 2.68 4,32 5697 7.62 10.91 14.20 17.49 20.78 
105.0 1,05 2.80 4,56 6.31 8.07 11.58 15.09 18.60 22,11 
110.0 1.06 2.94% 4,81 6.68 8.56 12.31 16.05 19.80 23,55 
115.0 1.08 3.08 5,08 7.09 9.09 13,09 17.10 21.10 25,12 
120.0 1.09 3.24 5,38 7.52 9.66 13.94% 18.23 22.51 26.79 
125.0 1.11 3.40 5.69 7.99 10.28 14.86 19.44 24,03 28.61 
130.0 1.13 3,58 6,03 8.49 10.94 15.85 20.76 25.67 30.58 


se aap ac cae a alee all 
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Table 2 (continued) 


WATER TEMPERATURE = 60.0 DEGREES FAHRENHEIT 


W=0 W=2.5 W=5.0 W=7.5 W=10 We15 W=20 W=25 W=30 
kq Temp 


40.0 88 1.78 2,67 3.57 4.47 6.26 8.05 9.85 11.64 
45.0 89 1,83 2.77 3.71 4.65 6.53 8.41 10.29 12,17 
50.0 91 1,89 2.88 3.87 4.85 6.82 8.80 10.77 12,74 
55.0 092 1.96 2,99 4.03 5.06 7.14% 9.21 11.28 13,35 
60.0 093 2.02 3,11 4.20 5.29 7.47 9.65 11.82 14,00 
65.0 095 2,09 3,24 4.38 5.53 7.82 10.11 12.40 14,70 
70.0 0960 2.17) = 3437 4.58 5.79 8.20 10.61 13,02 15,44 
75.0 298 2,25 3.52 4.79 6.06 8.60 11,14 13,69 16,23 
80,0 099 §=©2,.33 86367 5.01 6.35 9.03 11.71 14.39 17,07 
85.0 1.00 2.44% 3,87 5.30 6.74 9.60 12.47 15.33 18.20 
90.0 1.02 2.54 4,06 5.58 7.10 10.15 13.19 16,23 19,28 
95.0 1.03 2.65 4,27 5.89 7.56 10.74 13.97 17.21 20,44 
100.0 1.05 2.77 4,49 6.21 7.93 11.38 14.82 18.26 21.71 
105.0 1.06 2.90 4,73 6.57 8.40 12.07 15.74% 19.41 23,08 
110.0 1.08 3,04 4,99 6.95 8.91 12.82 126.7% 20.65 24.57 
115.0 1,09 3.18 5,27 72.36 9.45 13.63 17,81 21.99 26.17 
120.0 1.11 3.34 5,58 7.81 10,04 14.51 18.98 23.4% 27,91 
125.0 1.12 3.51 5,90 8,29 10.68 15.45 20.23 25.01 29.79 
130.0 1.14% 3.70 6,25 8.81 11.36 16.47 21.59 26.70 31.81 


WATER TEMPERATURE = 65.0 DEGREES FAHRENHEIT 


W=0 =W=2.5 W=5.0 W=7.5 W=10 W=15 W=20 W=25 W=30 
Eq Temp 


40.0 89 1.85 2,80 3.75 4.71 6,€1 8.52 10,42 12,33 
45.0 091 1.89 2.88 3.87 4.85 6,82 8,80 10.77 12,74 
50.0 092 1.96 2.99 4.03 5.06 7.14 9,21 11,28 13,35 
55.0 093 = 2,02 3.11 4.20 $2.29 7.47 9,65 11.82 14,00 
60.0 095 2.09 3.24 4.38 S.53 7,82 10.11 12.40 14,70 
65.0 296 2.17 3.37 4.58 5.79 38.20 10.61 13.02 15,44 
70.0 098 = 2425 3.52 4.79 6,06 8.60 11.14% 13.68 16,23 
75.0 099 2.33 3.67 5.01 6.38 9,03 11,71 14,39 17,07 
80.0 1.00 2.42 3.83 5.24 6.66 9.48 12,31 15.14 17,96 
85.0 1.02 2.51 4,00 5.49 6.98 9.97 12.95 15.93 18.92 
90.0 1.03 2.63 4,23 5.82 7.42 10.61 13.80 17.00 20.19 
95.0 1.05 2.748 4%, 44 6.13 7.83 11.22 14.61 18,01 21.40 
100.0 1,06 2.87 4,67 6.47 8.28 11.89 15.49 19,10 22.71 
105.0 1,08 3.00 4,92 6.84 @.76 12.60 16.44% 20.29 24,13 
110.0 1,09 3.14 5,19 7.23. 9.28 13,38 17,47 21.57 25.66 
115.0 1,11 3,29 5.48 7.66 9,85 14,21 18,58 22.95 27,32 
120.0 1.12 3.46 5.79 8.12 10.45 15.12 19.78 24,45 29,11 
125.0 1.14 3.63 6,12 8.61 11.11 16,09 21.07 26.06 31,04 
130.0 1,16 3.82 6.48 9.15 11.81 17.14 22.47 27,80 33,12 
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Table 2 (continued) 


WATER TEMPLRATURE = 70,0 DEGREES FAHRENHEIT 
W=0 W=2.5 W=5.0 W=7.5 W=10 W=15 W=20 W=25 W=30 


40.0 091 1.91 2.92 3.93 4.93 6.9% 8.95 10.97 12.98 
45.0 092 1.97 3.02 4.07 5.12 7.22 9.32 11.42 13.52 
50.0 093° 2.02 33411 4.20 5.29 7.47 9.65 11.82 14,00 
55.0 095 2.09 3,24 4.38 5.53 7,82 10.11 12.40 14.70 
60.0 096062417) 3,37 4.58 5.79 8.20 10,61 13.02 15.44 
65.0 098 =. 2.25 3.52 4.79 6,06 8,60 11.14% 13.68 16.23 
70.0 0990 24330 3.467 5.01 6.35 9,03 11.71 14.39 17.07 
75.0 1.00 2.42 3.83 5.24 6.66 9.48 12,31 15.14 17.96 
80.0 1.02 2.51 4.00 5.49 6.98 9.97 12.95 15.93 18.92 
85.0 1.03 2,61 4.18 5.76 7.33 10.48 13,63 16.78 19.93 
90.0 1.05 2.71 4.37 6.04 7.70 11,02 14.35 17.68 21.00 
95.0 1.06 2.84% 4,62 6.40 8.18 11.75 15,31 18.87 22.43 
1090.0 1.08 2.97 4,86 6.75 8.65 12.43 16.22 20.C0 23.79 
105.0 1.09 3.11 5.12 7.13 9.15 13,17 17.20 21.23 25.25 
110.0 1.11 3.25 5.40 7.54 9,68 13,97 18.26 22.55 26,84 
115.0 1.12 3.41 5.69 7.98 10.27 14,84 19,41 23.98 28.55 
120.6 1.14 3.58 6,01 8.45 10.89 15.77 20.6% 25.52 30.40 
125.0 1.15 3.76 6.36 8.96 11.57 16.77 21.98 27.18 32.39 
130.0 1.17 3.95 6.73 9.51 12.29 17.85 23.41 28.97 34.53 


WATER TEMPERATURE = 75.0 DEGREES FAHRENHEIT 


W=O0 W=2,5 W=5.0 W=7.5 W=10 W=15 W=20 W=25 wW=30 
Eq Temp 


40,0 092 1.99 3.05 4,11 5.18 
45.0 93 2.05 3.16 4,27 5.38 
50.0 095 2.11 3.27 4.43 =§,60 


7.3 9.43 11.56 13.69 

7.6 

7.9 
55.0 096 0=— 2,17 3.37 4.58 5.79 8.2 

8.6 

9.0 


1 

0 9.82 12.0% 14.27 
2 19,24 12.57 14.89 
0 10.61 13.02 15.44 
$0.0 098 82.25 3.52 4.79 6,06 O 11.14% 13.69 16.23 
65.0 099 2.33 3.67 5.01 6.35 3 11.71 14.39 17.07 
70.0 1,00 2042 3.83 5.24 6,66 9.48 12,31 15.14% 17.96 
75.0 1,02 2.91 4,00 $5.49 6.98 9.97 12,95 15.93 18.92 
80.0 1.03 2.61 4.18 5.76 7.33 10.48 13.63 16.78 19.93 
85.0 1,05 2071 4,37 6.04 7,70 11,02 14.35 17.68 21.00 
90.0 1.06 2.82 4.57 6.33 8.09 11.60 15.12 18.63 22.14 
95.0 1.08 2.93 4.79 6.65 8.50 12.2% 15.93 19.64 23.35 
100.0 1.09 3.0F 5.07 7.06 9,04 13,02 17.00 20.97 24.95 
105.0 1.11 3.22 5233 7.45 9,56 13,79 18.01 22.2% 26.47 
110.0 1,12 3.37 5.62 7.87 10,12 214,61 19.11 23.61 28,11 
115.0 1.14 3.53 5.93 8.32 10.72 15.51 20.30 25.08 29.87 
120.0 1.15 3.71 6.26 8.81 11,36 16.47 21.57 26.67 31.78 
125.0 1.17 3.89 6.61 9.33 12,06 17,50 22.°% 28.39 33.83 
130.0 1.19 4,09 6.99 9,90 12,80 18.61 24.42 30.23 236.04 
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Table 2 (continued) 


WATER TEMPERATURE = 80.0 DEGREES FAHRENHEIT 


WO = W=2.5 =W=5.0 W=7.5 W=10 W=15 W=20 W=25 W=30 
Lq Temp 


40.0 094% 2,06 3.19 4.32 5.45 7.70 9.96 1. °2 14,47 
4$.0 095 2.13 3,30 4.48 5.66 8.01 10.37 12,72 15,08 
50.0 96 8=62219 3.42 4.65 §.88 8.35 10.81 13,27 15.73 
55.0 098 = §=2426 3.55 4,84 6.13 8.71 11.29 13,86 16,44 
60,0 0990s 2.330 33,67 5.01 6.35 9.03 11.71 14,39 17,07 
65.0 1.00 2.42 3.83 5.24% 6.66 9.48 12,31 15,14 17,96 
70.0 1.02 2.51 4,00 5.49 6.98 9.97 12.95 15.93 18,92 
78.0 1.03 2.61 4,18 5.76 7,33 10.48 13.63 16.78 19,93 
80.0 1.05 2.71 4,37 6.04 7.70 11.02 14.35 17,68 21,00 
85.0 1.06 2.82 4.57 6.33 8.09 11.60 15.12 18,63 22.14 
90.0 1.08 2.93 4,79 6.64% 8,50 12.21 15.93 19,64 23,35 
95.0 1.09 3.05 5,01 6.98 8.94% 12.86 16.79 20.71 24,64 
100.0 1.11 3.18 5,25 7.33 9.40 13.55 17.70 21,84 25,99 
105.0 1.12 3.34 5,56 7.78 10.01 14.45 18.89 23,33 27.77 
110.0 1,14 3.50 5,86 8.22 10.58 15.30 20.02 24.75 29,47 
115.0 1,15 3.66 6,18 8.69 11.20 16.22 21.25 26.27 31.29 
120.0 1.17 3.84% 6,52 9.19 11.87 17.21 22.56 27.91 33.26 
125.0 1.18 4.03 6,88 9.73 12.58 18.28 23.98 29.68 35.37 
130.0 1.20 4.24 7,27 10.31 13.35 19.42 25.50 31.57 37.65 


WATER TEMPERATURE = 85.0 DEGREES FAHRENHEIT 


W=0 W=2,5 W=5.0 W=7.5 W=10 W=15 W=20 W=25 W=30 


Eq Temp 

40,0 095 2.15 3,35 4.54 5.74% 8.14% 10,53 12.93 15,33 
45.0 0960 26213 4G 4.71 5.96 8.46 10.96 13.46 15,96 
50.0 098 2.28 3.59 4.89 6.20 8.81 11.42 14,04 16,65 
55.0 0990 2.360 372 5.09 6.46 9.13 11.93 14,66 17,39 
60.0 1.00 2.44% 3,87 5.30 6.74 9.60 12.47 15.33 18,20 
65.0 1.02 2.51 4,00 5.49 6.98 9.97 12.95 15.93 18,92 
70.0 1.03 2.61 4,18 5.76 7.33 10.48 13.63 16.78 19,93 
75.0 1.05 2.71 4.37 6.0% 7,70 11.02 14.35 17.68 21,00 
80.0 1.06 2.82 4,57 6.33 8.09 11.60 15,12 18.63 22,14 
85.0 1.08 2.93 4.79 6.64 8,50 12.21 15.93 19,64 23,35 
90.0 1.09 3.05 5,01 6.98 8.9% 12.86 16,79 20,71 24.63 
95.0 4.12 3.18 5,25 7.33 9.40 13.55 17.70 21.84 25,99 
100.0 1.12 3.31 5.51 7.70 9,89 14.27 18.66 23.04 27,43 


105.0 1,14 3.45 $277 8.09 10.41 15.04 19,68 24,31 28,95 4 
110.0 1.15 32.63 6,11 8.60 11.08 16.04% 21,00 25.97 30.93 ; 
115.0 1,17 3.81 &,44 9,08 11.72 16.99 22,27 27.54% 32,82 
120.0 1.18 3299 6.79 9.60 12.40 18.01 23.62 29.24 34,85 

125.0 1.20 4.19 7.17 10.16 13.14% 19.11 25,08 31.05 37,03 
130.0 1.22 4.40 7.57 10.75 13.93 20.29 26,65 33,01 39,37 ai 
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Table 2 (continued) 


WATER TEMPERATURE = 90.0 DEGREES FAHRENHEIT 


W=0 W=2.5 W=5.0 W=7.5 W=10 W=15 W=20 W=25 


Eq Temp 
40.0 96 2.24 3.51 4.79 «©6606 8.61 11.16 13,71 
45.0 098 2.31 3.64 4.96 6.29 8.95 11.61 14,27 
50.0 99° 2.38 8=—33777 £015 6.54% 9.32 12,10 14,87 
55.0 1.00 2.46 3.91 5.36 6.f1 9.72 12.62 15,53 
60.0 1.02 2.54% 4,06 5.58 7.10 10.15 13.19 16,23 
65.0 1.03 2.63 4,23 5.82 7.42 10.61 13.80 17,00 
70.0 1.05 2.71 4,37 6.0% 7,70 11.02 14.35 17.68 
75.0 1.06 2.82 4,57 6.33 8.99 11.60 15.12 18.63 
80.0 1.08 2.93 4.79 6.64 8,50 12.21 15.93 19,64 
85.0 1.09 3.05 5.01 6.98 8.94% 12.86 16.79 20.71 
90.0 1.11 3.18 5.25 7.33 9.40 13.55 17.70 21,84 
95.0 1.12 3.31 5.51 7.70 9.89 14.27 18.66 23,04 
100.0 «62.44% 3.45 5.77 8.09 10.41 15.04 19.68 24,31 
105.0 1.15 3.60 6,05 8.50 10.95 15.85 20.75 25,65 
110.0 1.17 3.76 6,35 8.94% 11.53 16.71 21.89 27,07 
115.0 1.18 3.96 6.73 9.50 12.27 17.82 23.36 28.91 
120.0 1.20 4.14 7.09 10.04% 12,98 18.87 24.76 30.65 
125.0 1.22 4.35 7.48 10.61 13,74 20.00 26.27 32,53 
130.0 1.23 4.56 7.89 11.23 14,56 21.22 27.88 34,54 


W=30 


16.26 
16.93 
17.65 
18,43 
19,28 
20,19 
21,00 
22014 
23.35 
24.63 
25.99 
27.43 
28,95 
30.55 
32.25 
34.45 
36.55 
38.79 
41.21 
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Appendix A 


TEMPERATURE DEPENDENCE OF RELATIVE ISOBARIC SLOPES 


To include relative isobaric slopes in the vertically averaged 


momentum equations, replace the terms aes and . by the expanded 


terms g(dE/dax + < i > WD and g(d€/dy + < i > 7 respectively. 


XopsE YoPs& 
Pritchard gives the following expression for the vertically averaged 
; F (19) 
isobaric slope: 
‘ 1 3p H 30 as 
= W/2H-S = SCT 
nr) h / ao ox 2p 96 dx (Al) 


where 0 = vertically averaged density, and s = vertically averaged 
salinity. In terms of the specific volume defined by a = 1/0, this 


becomes 


i, = (H/2) @s/3x) (- + 32, : (A2) 


The quantity (- + 3% is the coefficient of saline contraction and 


may be obtained according to Eckart from the Tumlirz equation of 


state given by (29) 


(p + P,) (a - a) = (A3) 


where T is in °C and 


1779.5 + 11.25T - .0745T" - (3.80 + .01T)s , 


> 
tt 


a = .6980 , 
fo) 
P= 5890 + 38T - 375 T? + 3s , 
P = total pressure (in atmospheres) , 


it 


a = specific volume in ml/gm . 


For H < 100 m we can ignore p compared to Por For example, at 10 
meters depth, P/P, is less than to: Since the empirical formula 
is accurate only to about 1%, our approximation is well justified. 
Evaluation of the saline contraction is now straightforward and 


the result is: 


a a 
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3.8 + .01T + 3 (\/p_) 
i 5. (2) ‘C8 fx) (A4) 


x,h A+ a Po 
A similar result applies for x replaced by y. 
If data is given in terms of chlorinity, one must use the fol- 
lowing equation relating salinity and chlorinity (both in parts per 


thousand) ; 617) 


s = .03 + 1.8050 Cl . (A5) 


Chlorinity is defined as the chlorine equivalent in parts per thousand 
of all precipitated halides. Equation (A5) does not apply to estuaries 


whose constituents are altered by significant fresh water flow. 


Let s' and Cl’ denote the salinity and chlorinity, respectively, 


es pressed in g/liter sea water. Then 


Ss =o0s ’ 
(A6) 
Cl' = pCl , 
and 
s' = .030 + 1.805 Cl’ 
(A7) 
= ees ee + 1.805 cl' 
a + X;p 
° fo) 
In terms of s', i becomes 
x,h 
i, = (H/2) (as"/ax) (- +2%,) . (A8) 
x,h ; ads! 
We also obtain 
Se Oe es oe ert, Se 1 
~ @ 9s! a(as'/da) ds s s ‘ (A9) 
—-= — +s 
aa a aa 
and 
as'/dx = .03(3p/dx) + 1.805(sC1'/5x) . (Al0a) 


5 
In (Al0a) the first term on the right hand side is roughly 10 


times snaller than the second term, so we can write approximately 
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ds'/dx = 1.805(9C1'/dx) (A10b) 


In terms of chlorinity, the isobaric slope is given by 


i = 1.805(H/2) {(aci' /3x)/ 
x,h 
P 
-03 () 

1.805 cl1' + 4 

| A + a, 3.8 + ,O1T + EL | (All) 
Po Py 
where s appearing in the definitions of i and Po is given by 
cl' ’ 

s = .03 + 1.805 T.o- -03 + 1.770 Cl (A12) 


In (Al2) we have used the value p = 1.02, which is correct within 
1% for the range of temperature and salinity under consideration. 
Errors of order 1% in s will produce errors of order only .001% in 
A and rs Similar arguments lead us to drop the term -03(a, + Vi ee 
in (All). 

Temperature gradients will also contribute to the relative iso- 


baric slope. For example, (Al) is modified according to 


: 19p as, lap at 
th tay) (#22 ax * paT ax} * 2) 


The coefficient of thermal expansion is not given as accurately by 
our equation of state as is the saline contraction. However, this 
is compensated by the smaller magnitude of the former effect. 


From (A3) we find 


oer) eee) ee Os a 
0 aT a aT a [dT Arp.) 
Q/p.) €p_/3T) - (8d/aT) 
Shp fe Oe ON (A14) 
A+ ap 
oe) 
Q/p)) (38 - .750T) - (11.25 - .149T - .018 C1") 
A + oP. 


Finally, the relative isobaric slope as a function of chlorinity 
{in g/liter), temperature (°C), the corresponding gradients, and 


depth is given by 
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; 1d acl’ , Lap ar 
do PANS OGlt Oe oo aT =| (Al5) 
: ' 
= 1.805(H/2) a _C1'/dx - (A16) 
1.805Cl' + 2 


3.8 + .01T + (3\/p.) 


+ (H/2) (dT/dx) 


Q/p,) (38 - .750T) - (11.25 - .149T - .918 Cl") 
A + oP 


where A, Po? and a, are defined in (A3) and the value of s used to 
evaluate A and Py is given by (Al12) 

Equation (Al6) may be used for machine computation when large 
salinity or temperature gradients make it desirable to employ local 
values of the saline contraction or thermal expansion coefficient. 
In most cases it will be better to use averaged values for the 
coefficients over the entire field. Typical numerical results for 
thermal and chlorinity coefficients are 2x10°4 and 1.3x10"7 respec~ 


tively. 
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Table A 


Table of Z(T,C) where Z is the coefficient of chlorinity 
contraction as a function of centigrade temperature 
and chlorinity in grams per liter. 


Chlorinity in g/l 
10,000 
10.250 
10,500 
10.750 
11.000 
11.250 
11.500 
11.750 
12,000 
12.250 
12.500 
12.750 
13,000 
13,250 
13.500 
13.750 
14.000 
14.250 
14,500 
14.750 
15,000 
15.250 
15.5006 
15.750 
16.000 
16.250 
16.5900 
16.750 
17,000 
17,250 
17,590 
17.750 
18.000 
18.250 
18.500 
18.750 
19.000 
19.250 
19.500 
19.750 
20,000 


1.350-03 
1,349-03 
1.348-03 
1.347-03 
1. 346-03 
1.345-03 
1. 344-03 
1.343-03 
1.342-03 
1.341-03 
1,340-03 
1.339-03 
1.338-03 
1.337-03 
1.336-03 
1.335-03 
1.334-03 
1.333-03 
1.332-03 
1.331-03 
1.330-03 
1.329-03 
1,.328-03 
1.327-03 
1.326-03 
1.325-03 
1.324-03 
41.323-03 
1.322-03 
1.321-03 
1.320-03 
1,319-03 
1.318-03 
1.317-03 
1,316-03 
1,315-03 
1.315-03 
1,314-03 
1.313-03 
1,.312-03 
1.311-03 


T=15 


1.335-03 
1.334-03 
1.333-03 
1.332-03 
1,.331-03 
1.330-03 
1.329-03 
1.328-03 
1.327-03 
1.326-03 
1.325-03 
1,324-93 
1,323-03 
1.322-03 
1.321-03 
1.320-03 
1.319-03 
1,318-03 
1.317-03 
1.316-03 
1,315-03 
1.314-03 
1.313-03 
1.312-03 
1,312-03 
1,311-03 
1,310-03 
1, 309-03 
1.308-03 
1,307-03 
1.306-03 
1,305-03 
1,304-03 
1,303-03 
1,302-03 
1, 301-03 
1,300-03 
1,299-03 
1,298-03 
1,297-03 
1,.296-03 


T=20 


1,324-03 
1.323-03 
1,.322-03 
1,.321-03 
1.320-03 
1,319-03 
1.318-03 
1.317-03 
1.316-03 
1.315-03 
1,.314-03 
1.313-03 
1,313-03 
1.312-03 
1.311-03 
1.310-03 
1.309-03 
1.308-03 
1.307-03 
1.306-03 
1.305-03 
1.304-03 
1.303-03 
1,302-03 
1.301-03 
1. 300-03 
1.299-03 
1,.298-03 
1.297-03 
1.297-03 
1.296-03 
1,295-03 
1, 294-03 
1.293-03 
1.292-03 
1,291-03 
1,290-03 
1,.289-03 
1,288-C3 
1.287-03 
1.286-03 


1,317-03 
1.317-03 
1.316-03 
1.315-03 
1.314-03 
1.313-03 
1,312-03 
1.311-03 
1.310-03 
1.309-03 
1.308-03 
1.307-03 
1.306-03 
1.305-03 
1.304-03 
1.303-03 
1.302-03 
1.301-03 
1.301-03 
1.300-03 
1.299-03 
1.298-03 
1.297-03 
1.296-03 
1,295-03 
1,2¢4-03 
1.283-03 
1.292-03 
1,291-03 
1,290-C3 
1.289-03 
1,289-03 
1,288-03 
1,287-03 
1.286-03 
1.285-03 
1, 284-03 
1,283-03 
1,282-03 
1.281-03 
1.280-03 


T=30 


1,315-03 
1,314-03 
1,313-03 
1,312-03 
1,.311-03 
1,310-03 
1,309-03 
1,368-03 
1.307-03 
1.306-03 
1.305-03 
1,.305-03 
1, 304-03 
1,303-03 
1,302-03 
1,301-03 
1.300-03 
1.299-03 
1,298-03 
1,297-03 
1,296-03 
1,295-03 
1, 294-03 
1,294-03 
1,293-03 
1,292-03 
1,291-03 
1,290-03 
1,289-03 
1.288-03 
1,287-03 
1,286-03 
1,285-03 
1,284-03 
1,284-03 
1,283-03 
1,282-03 
1,281-03 
1,280-03 
1,279-03 
1, 278-03 


a oeeeneeeenmnaneenennmnnamnennenemmnemnenimn tas eye eee nes Ae 
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Appendix B 


MODIFIED DEWPOINT TEMPERATURE 


Rather than regard Te as a function of Q = absorbed radiation, 


W, T. and T, we can take advantage of the following procedure which 


d 
considerably simplifies interpolation in the tables. The air tem- 
perature and dewpoint temperature enter the equation for a5 only in 


the term (e - e + BT_ - BT_) where, as before, e = e(T,). Let 
eq a eq a a d 


us define a "modified dewpoint cenpereeuretes J, according to 
- ~- = - + pe - 
Ces e, + BCT. T)) San e(J) BC n J) (Bl) 
or, 
e(J) = BJ = e, + BT. (B2) 


The modified dewpoint J may be found by Newton's method and tabulated 


as a function of T and T The equilibrium temperature may then 


dad’ 
be found as a function of J, W and Q. In particular, 


T (T.> Ty W, Q) 


> » W, B 
- ES Gage) (B3) 


e 


so that we need only the first section (where V = 0 and T = R) of 
Table 3. 


J is tabulated in Table B as a function of tT. and V = t. - T,. 
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Table B 


Table of modified dewpoint temperature J(T,R) as a function 
of air temperature T and v= [air temperature - dewpoint 
temperature], all in degrees Fahrenheit, 


v=0 V=5 V=10 V=1S V=20 V=25  v=30  v=Ho 


Ts 

QO. 29 ~.8 “1.4 1.9 -2.3 -2.6 -2,8 -3.2 
5.6 5.1 4,2 3.4 2.8 2.3 1.9 1.6 1.2 
10. 10.1 8.0 6.1 7.4 6.8 6.3 5.9 5.4 
15. 15.1 13.9 12,8 11.9 11,2 10.7 10,2 9.5 
20. 20.1 18.7 17.5 16.5 15.6 14.9 14.4 13.6 
25. 25.1 23.5 22.1 20.9 20.0 19.2 18.5 17.5 
30. 30.2 28.3 26.7 25.4 24.2 23.3 22,6 21,4 
35. 35.2 33.1 31,3 29.8 28.5 27.4 26.5 25.2 
40, 40.2 37.9 35.9 34.2 32.7 31.5 30.5 28.9 
45, 45.2 42.7 40.5 38.6 37.0 35.6 34,4 32.6 
59, 50.0 47.5 45.1 43.0 41.2 
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